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We show that the fundamental seismic shear mode, observed as a quasi-periodic osciUation in 
giant flares emitted by highly-magnetized neutron stars, is particularly sensitive to the nuclear 
physics of the crust. The identification of an oscillation at « 30 Hz as the fundamental crustal shear 
mode requires a nuclear symmetry energy that depends very weakly on density near saturation. If 
the nuclear symmetry energy varies more strongly with density, then lower frequency oscillations, 
previously identified as torsional Alfven modes of the fluid core, could instead be associated with 
the crust. If this is the case, then future observations of giant flares should detect oscillations at 
around 18 Hz. An accurate measurement of the neutron skin thickness of lead will also constrain 
the frequencies predicted by the model. 

PACS numbers: 26.60.Gj,97.10.Sj,26.60.-c,21.65.Ef 



Giant x-ray flares from highly magnetized neutron 
stars are powered by catastrophic reconfigurations of the 
decaying field. Pinning of the rapidly evolving field to the 
solid crust triggers an associated starquake and generates 
global seismic vibrations [l[ , detectable as quasi-periodic 
oscillations (QPOs) in the x-ray afterglow [1, 0, IjiSl- 

The current picture is that the QPOs, which range in 
frequency from 18 Hz up to 1800 Hz, result from tor- 
sional (twisting) motions of the star. The lowest (18, 26 
Hz) frequencies were initially associated with the mag- 
netized fiuid core, while the other (> 28 Hz) frequencies 
were interpreted as shear modes of the solid crust. The 
problem is actually more complicated, as the fluid core 
admits an Alfven continuum, and crust and core arc cou- 
pled by the strong magnetic field. However even in this 
more complex system frequencies close to the pure crustal 
frequencies can still emerge, from global modes in which 
crust motion dominates [y, 0, @1 ■ If this picture is cor- 
rect, then the observed flare QPO frequencies probe the 
shear properties of the neutron star crust. This is partic- 
ularly exciting because the frequencies can be measured 
relatively accurately, to within a few percent or better. 

In this paper, we calculate the frequency of shear os- 
cillations of the neutron star crust, and show that they 
depend sensitively on a particular aspect of the nuclear 
physics input: the nuclear symmetry energy. The sym- 
metry energy, the energy cost of creating an isospin asym- 
metry in nucleonic matter, is one of the most si gnifi cant 
uncertainties in the description of the crust [ol. lid. [lH . 
The sensitivity of the flare QPO frequency to the symme- 
try energy implies that they can constrain the properties 
of the nuclear symmetry energy. We also make a novel 
connection to the proposed measurement of the neutron- 
skin thickness of a lead nucleus (the difference between 
the neutron and proton radii), to be performed next year 



at Jefferson Lab in the PREX [13, [ist experiment. 

At lower densities, the outer crust of the neutron star 
consists of nuclei embedded in a sea of degenerate elec- 
trons. As density rises, nuclei become progressively heav- 
ier and more neutron-rich. Once density increases past 
the neutron drip point, 4 x 10^^ g/cm^, it becomes en- 
ergetically favourable for neutrons to drip out of nuclei. 
This begins the inner crust region of the neutron star, 
where exotic nuclei are embedded in a sea of supcrfluid 
neutrons. At a density of around 1.5 x 10^^ g/cm'^, nuclei 
arc no longer favored and dissolve into their constituents. 

For the model of the inner neutron star crust, wc em- 
ploy the liquid droplet model from Ref. llj, as updated 
in Ref. [3l • This model consistently describes the nuclei 
and the dripped neutrons in the crust with one equation 
of state for homogeneous nucleonic matter, and is very 
similar to that employed for the description of matter in 
core-collapse supernovae [l^. For the homogeneous nu- 
cleonic matter EOS, we use the Skyrme model which 
enables us to describe matter over a large enough density 
range to obtain masses and radii, and also can describe 
neutron matter at very low densities in the inner crust. 

The determination of the crust-core transition density 
is difficult because the relevant energy surfaces are very 
flat near the transition, so we choose to fix the transition 
density at 0.07 fm~'^. Increasing this value to 0.1 fm~'^ 
decreases the frequency of the radial overtones (which de- 
pend strongly on crust thickness), but changes the funda- 
mental crust mode frequency only by a percent or so. We 
choose to use the outer crust model from Ref. [l7|, but 



have checked that using alternate models from Ref. [18| 
does not change the inferred QPO frequencies (the shear 
mode frequencies are largely insensitive to uncertainties 
in the neutron drip point [19|). 
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The shear modulus of the crust [20| is 
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Here Z is the atomic number of the ions, rii is the ion den- 
sity, and a = {S/inni)^^^ is the average inter-ion spacing. 
The parameter F = {Ze)^ / aksT , where T is the temper- 
ature, measures the ratio of the Coulomb to thermal en- 
ergies, and we use F = Fn = 173 to determine the upper 
boundary of the crust [21|. Equation ([1]), derived using 
a Monte Carlo simulation of the Coulomb interactions in 
a neutron star crust, assumes that the contribution to 
the shear modulus from neutron-lattice interactions and 
neutron-neutron interactions is small. 

Neutron stars and their magnetospheres admit many 
types of vibration, driven by different restoring forces. 
The identification of the QPOs excited in magnetar flares 
with crust shear modes is based on several factors. Shear 
modes have a lower excitation energy (compared to vi- 
brations involving bulk compression), and damp suffi- 
ciently slowly to explain the observed QPO durations 
[2^ . Excitation of crust motion is particularly likely if, 
as argued by [l^l , the yielding of a strained crust triggers 
the flares. Coupling of a twisting crust to the external 
field also provides a plausible mechanism for modulating 
the x-ray emission [3|. Most importantly the observed 
QPOs match expectations from models in terms of both 
frequency and the scaling between different harmonics. 
Coupling of crust and core by the magnetic field does 
complicate the system but, as demonstrated by Q, does 
not prevent the emergence of the natural (uncoupled) 
frequencies. We can therefore compute crust shear mode 
frequencies assuming free slip between crust and core. 

We employ the simple Newtonian perturbation model 
for torsional shear oscillations used by [l^, H^. This 
model uses a plane-parallel geometry with constant grav- 
itational acceleration, computed beforehand using the 
TOV equations. The Newtonian equations of hydrostatic 
equilibrium then determine the crust density profile. Us- 
ing a slab rather than spherical geometry allows us to 
incorporate the magnetic field without difficulty; we as- 
sume a constant field B = Bz. In order to correct for 
the error that this introduces, and to ensure that we re- 
cover the correct spheric al g eometry frequencies in the 
zero field limit we follow [26| and mimic a spherical ge- 
ometry by setting V^^^ = -[(/ + 2){l- l)/i?^]^, I being 
the standard angular quantum number. 

For pure toroidal shear modes, which are incompress- 
ible and have no vertical component of displacement, the 
perturbation equations for the horizontal displacement ^ 
reduce to 
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FIG. 1: Shear speeds Vs (solid lines) and the neutron skin 
thickness of lead 5R (dashed lines) as a function of the sym- 
metry energy. The abcissa is the magnitude of the symmetry 
energy at the nuclear saturation density and the ordinate is 
the derivative of the symmetry energy. 



where va = B /{inp)^^^ is the Alfven speed. The shear 
speed is Vs = (ii/p)^^^. We assume a periodic time de- 
pendence exp(iwt), uj being the frequency, and correct 
for gravitational redshift to obtain observed frequency. 
Primes indicate derivatives with respect to z. 

Using the schematic model of Ref. [ll[, we can con- 
struct a neutron star crust from an arbitrarily chosen 
symmetry energy. This crust model includes some cor- 
rections to the nuclear masses due to the surrounding 
neutron gas. Fig. [T] displays the connection to the nu- 
clear symmetry energy as a function of density, i5sym('^), 
relating the magnitude of the symmetry energy, S, at the 
nuclear saturation density, no = 0.16 fm~"^, and a param- 
eter related to the derivative of the symmetry energy, 
L = 'SnQdEsym/dn\n=nai to the shear speed, Vs- The 
shear speed is large if the magnitude of the symmetry 
energy is large and the density dependence is small. Also 
displayed is the relationship to the neutron skin thick- 
ness of lead, ((5 i?), from the Typel-Brown correlation 27|, 
which depends most strongly on L. 

For a full oscillation model, we need a complete model 
for hadronic matter up to the densities in the center 
of neutron stars. The shear speeds as a function of 
mass density are given in Fig. [2] for the Skyrme mod- 
els, BSkl4 ji^, Gs [2^, NRAPR [l], Rs [2^, SLy4 
SkI6 SkOpi, SkT2 We also use the model 
from APR [3J|, for which SR has been estimated at 
0.20 fm loll. Previous calculations such as those in 
Refs. [24, bdl were based on the crust of Ref. [351. This 
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FIG. 2; The shear speeds as a function of density for the crust 
models used in this work, and the shear speed of Ref. [35| 
(labeUed DHOl), which was based on SLy4. Also shown is 
the neutron skin thickness of lead, 5R, for each model. 
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FIG. 3: The crust oscillation frequencies as a function of neu- 
tron star mass, for both the fundamental (n = 0, Z = 2) tor- 
sional shear mode and the first radial (n = 1) overtone. The 
curves end at the maximum mass. The arrows on the right in- 
dicate QPO frequencies measured during the 2004 hyperflare 
from SGR 1806-20 UM^- 



also used the Skyrme model SLy4 [30[ and is shown as a 
dotted Une, labelled DHOT The difference between that 
result and our model, derived from the same Skyrme in- 
teraction, provides an estimate of the uncertainties in the 
construction of the masses of the neutron-rich nuclei in 
the crust, given the same underlying nucleon-nucleon in- 
teraction. That difference is smaller than the variation 
arising from the nucleon-nucleon interaction itself, which 
is as much as a factor of 4. This large variation in the 
shear speed is the source of the sensitivity of the QPO 
frequencies to the nuclear physics input. As expected 
from Fig. [H there is a correlation between the neutron 
skin thickness in lead, with smaller skins generally giv- 
ing larger shear speeds. The upcoming PREX measure- 
ment will help constrain the frequencies predicted by the 
model, reducing the uncertainty originating from our lack 
of knowledge of the symmetry energy. 

Figure [3] shows the frequencies of the fundamental 
n = 0, Z = 2 crust shear mode and the first (n = 1) 
radial overtone as a function of neutron star mass for 
the various crust models. Also shown in this figure are 
some of the QPO frequencies measured during the 2004 
hyperflare from SGR 1806-20. The dramatic effect of un- 
certainty in the symmetry energy on shear speed is evi- 
dent in the spread of mode frequencies. The frequency of 
the fundamental mode is larger for smaller neutron skin 
thicknesses, however all of the models have difficulty ex- 
plaining a fundamental above 22 Hz. Previous studies 
including which interpreted the 28 Hz QPO in SGR 
1900+14 and the 29 Hz QPO in SGR 1806-20 as the 
fundamental, employed the higher shear speed model of 
[sit , thereby yielding higher frequencies. Once the cor- 



rect I scaling is taken into account, however, even this 
crust model would necessitate a rather low stellar mass 
[2^. The models that we have considered, which have 
a symmetry energy that varies more strongly with den- 
sity, cannot explain a fundamental frequency at 28-29 
Hz. The only observed QPO that would fall into the 
predicted range is the 18 Hz QPO detected in the SGR 
1806-20 hyperflare 0, Q and previously interpreted as a 
torsional mode of the magnetized fluid core 

Rclativistic perturbation calculations predict very sim- 
ilar frequencies for the fundamental to those generated 
by Newtonian calculations [2^. The mode frequencies 
are also largely insensitive to magnetic field effects, since 
the electron Fermi momentum in the inner crust is much 
larger than the energy spacing between Landau levels 
except at fields larger than the 10^"'" ^'^ G implied for 
magnetars. Note also that the fundamental frequency is 
nearly mass- independent, which means that uncertain- 
ties in the neutron star mass will not spoil the connec- 
tion between the frequency and the nuclear symmetry 
energy. It is not well known how finite-size effects and 
the presence of extremely-deformed nuclei (i.e. nuclear 
pasta) may affect the shear modulus. Recent calculations 
by [36| , which re-examine Equation ([T]) would exacerbate 
the problem since they indicate that shear modulus is 
actually lower by ^ 10%. This would reduce mode fre- 
quencies by a further 3-4%. Finally, the structure of the 
neutron-rich nuclei are not well understood and this also 
adds a systematic uncertainty to our results. 

Our crust model thus suggests that a revision of the 
previous interpretation of the 28 Hz mode as the funda- 
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mental shear mode may be necessary. One alternative is 
that the fundamental mode may have a lower frequency 
(an idea also explored in [s^])- For this to be plausible, it 
should be possible to fit the sequence of QPOs detected 
in each flare with the expected scaling in L In the SGR 
f806-20 giant flare, a QPO was detected at f8 Hz. If 
this is the fundamental (/ = 2), then a reasonable fit to 
the higher frequencies can be obtained with the follow- 
ing mode identifications: 29 Hz (/ = 3), 92 Hz (/ = 10), 
150 Hz (/ = 16). The one QPO that would not fit the 
sequence is the 26 Hz QPO, but this could also be accom- 
modated if the fundamental were even lower, at « 11 Hz. 
Such low values are predicted by some of our models, al- 
though it might then be difficult to explain the 625 Hz 
QPO as an n = 1 radial overtone. Alternatively, the 26 
Hz QPO is a core-dominated mode. 

No lower frequency QPOs were detected in the SGR 
19004-14 hyperflarc. However the sequence of detected 
QPOs would be compatible with a fundamental at « 18 
Hz with the following mode identifications: 28 Hz (/ = 3), 
53.5 Hz {I = 6), 84 Hz {I = 9) and 155 Hz (/ ^ 17). As 
for SGR 1806-20, a lower frequency fundamental at w 11 
Hz could also be accommodated. With this in mind, 
we performed a thorough rotational phase and energy- 
dependent search of the data from the Rossi X ray Tim- 
ing Explorer for the SGR 1900+14 hyperflare to search 
for lower frequency QPOs. No detections were made, but 
due to poorer data quality we could not search for QPOs 
as weak as those detected in SGR 1806-20. At 18 Hz, 
for example, we were able to set a 3 sigma upper limit 
on a QPO amplitude of 7% rms: in the SGR 1806-20 
hyperflare the 18 Hz QPO had an amplitude of 4%. 

It is interesting that an 18 (or 11) Hz fundamental fre- 
quency might be able to explain both mode sequences. If 
this is a common feature then we would expect to detect 
such a feature in future giant flare lightcurves. The low- 
ering of crust shear mode frequencies and their proximity 
to the expected frequencies of magnetic torsional modes 
does however complicate efforts to use magnetar QPOs 
to measure interior field strengths 0, [s^l- A sequence 
of modes, including overtones, is likely to be necessary 
to enable correct mode identification. There are several 
other possible effects on the shear modulus, including 
that of shell effects, frozen impurities 38, 3§|, superfluid- 
ity and entrainment (4o| , and the magnetic field which we 
have not addressed. A more definitive mode assignment 
will have to wait until these effects are better understood, 
or more data is obtained. 

AWS was supported by JINA at MSU under NSF- 
PFC grant PHY 08-22648, and by NASA ATFP grant 
NNX08AG76G. ALW thanks JINA for hospitality. We 
also thank Ed Brown, Lars Samuelsson and participants 
of the University of Washington Institute for Nuclear 
Theory workshop "The Neutron Star Crust and Surface: 
Observations and Models" for stimulating discussions. 



[1 
[2 
[3 

[4; 

[s: 
[6: 

[7: 

[9: 

[lo; 

[11 

[12 
[13 
[14 



[is: 

[16: 
[17: 

[is: 

[19 



[20 

[21 

[22' 
[23 

[24 
[25 
[26 

[27 

[28: 

[29 

[30 
[31 

[32: 

[33 
[34: 



R. C. Duncan, Astrophys. J. Lett. 498, 45 (1998). 

G. L. Israel et al, ApJ 628, L53 (2005). 

T. E. Strohmayer and A. L. Watts, ApJ 632, Llll 

(2005) . 

A. L. Watts and T. E. Strohmayer, ApJ 637, LI 17 

(2006) . 

T. E. Strohmayer and A. L. Watts, ApJ 653, 593 (2006). 

K. Glampedakis, L. Samuelsson, and N. Andersson, MN- 

RAS 371, L74 (2006). 

U. Lee, MNRAS 374, 1015 (2007). 

Y. Levin, MNRAS 377, 159 (2007). 

A. W. Steiner, M. Prakash, J. M. Lattimer, and P. J. 
Ellis, Phys. Rep. 411, 325 (2005). 

K. Oyamatsu and K. lida, Phys. Rev. C 75, 015801 

(2007) . 

A. W. Steiner, Phys. Rev. C 77, 35805 (2008). 

C. J. Horowitz, S. J. Pollock, P. A. Souder, and 

R. Michaels, Phys. Rev. C 63, 025501 (2001). 

R. Michaels, P. A. Souder, and G. M. Urcioh (2005), 

http: / /hallaweb.jlab.org/parity/prex. 

S. R. Souza, A. W. Steiner, W. G. Lynch, 

R. Donangelo, and M. A. Famiano, (available at 

http://arXiv.org/abs/0810.0963). 

J. M. Lattimer and F. D. Swesty, Nucl. Phys. A 535, 331 
(1991). 

T. H. R. Skyrme, Nucl. Phys. 9, 615 (1959). 

G. Baym, C. Pethick, and P. Sutherland, Astrophys. J. 

170, 299 (1971). 

S. B. Riister, M. Hempel, and J. Schaffncr-Bielich, Phys. 
Rev. C 73, 35804 (2006). 

See EPAPS Document No. [/number will be in- 
serted by publisher/] for variations in the transi- 
tion density. For more information on EPAPS, see 
http://www.aip.org/pubservs/epaps.html. 
T. Strohmayer, H. M. van Horn, S. Ogata, H. lyetomi, 
and S. Ichimaru, ApJ 375, 679 (1991). 
R. T. Farouki and S. Hamaguchi, Phys. Rev. E 47, 4330 
(1993). 

R. C. Duncan, ApJ 498, L45 (1998). 

C. Thompson and R. C. Duncan, MNRAS 275, 255 

(1995). 

A. L. Piro, ApJ 634, L153 (2005). 

A. L. Watts and S. Reddy, MNRAS 379, L63 (2007). 

L. Samuelsson and N. Andersson, MNRAS 374, 256 

(2007). 

S. Typel and B. A. Brown, Phys. Rev. C 64, 027302 
(200i). 

S. Goriely, M. Samyn, and J. M. Pearson, Phys. Rev. C 
75, 064312 (2007). 

J. Friedrich and P.-G. Reinhard, Phys. Rev. C 33, 335 
(1986). 

E. Chabanat, Ph. D. Thesis, University of Lyon (1995). 
P.-G. Reinhard and H. Flocard, Nucl. Phys. A 584, 467 
(1995). 

P.-G. Reinhard, D. J. Dean, W. Nazarewicz, 
J. Dobaczewski, J. A. Maruhn, and M. R. Strayer, 
Phys. Rev. C 60, 014316 (1999). 

F. Tondeur, M. Brack, M. Farine, and J. M. Pearson, 
Nucl. Phys. A 420, 297 (1984). 

A. Akmal, V. R. Pandharipande, and D. G. Ravenhall, 
Phys. Rev. C 58, 1804 (1998). 



5 



[35] F. Douchin and P. Haensel, A&A 380, 151 (2001). 

[36] C. J. Horowitz and J. Hughto, (available at 

littp://arxiv.org/abs/0812.2650). 
[37] H. Sotani, K. D. Kokkotas, and N. Stergioulas, MNRAS 

375, 261 (2007). 



[38] P. B. Jones, MNRAS 306, 327 (1999). 
[39] P. B. Jones, MNRAS 321, 167 (2001). 
[40] N. Andersson, K. Glampedakis, and L. Samuelsson, MN- 
RAS 396, 894 (2009). 



